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Abstract. We review our studies on the p and a resonances properties derived from the Inverse Amplitude Method. In 
particular, we study the leading \/N c behavior of the resonances masses and widths and their evolution with changing m n . 
The 1 /N c expansion gives a clear definition of qq states, which is neatly satisfied by the p but not by the cr, showing that 
its dominant component is not qq. The m K dependence of the resonance properties is relevant to connect with lattice studies. 
We show that our predictions compare well with some lattice results and we find that the pniz coupling constant is m n 
independent, in contrast with the anil coupling, that shows a strong m n dependence. 
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Light hadron spectroscopy lies beyond the realm of perturbative QCD. At low energies, however, one can use the 
QCD low energy effective theory, named Chiral Perturbation Theory (ChPT) [1], to describe the dynamics of the 
lightest mesons. ChPT describes the interactions of the pseudo Nambu-Goldstone bosons of the QCD chiral symmetry 
breaking, namely, the pions, by means of a effective lagrangian compatible with all QCD symmetries involving only 
the pion field. The infinite tower of terms in this lagrangian is organized as a low energy expansion in powers of p 2 /A^, 
where p stands either for derivatives, momenta or masses, and A^ ~ 4-nf K , where f K denotes the pion decay constant. 
ChPT is renormalized order by order by absorbing loop divergences in the renormalization of higher order parameters, 
known as low energy constants (LECs), that parametrize the high energy QCD dynamics and carry no energy or mass 
dependence. They depend on a regularization scale fi but after renormalization the observables are independent of this 
scale. The value of the LECs depend on the underlying QCD dynamics and are determined from experiment. Up to 
the desired order, the ChPT expansion provides a systematic and model independent description of how observables 
depend on some QCD parameters like the light quark mass m = (m u + m<j/2) or the number of colors, A^ c [2] . 

The use of ChPT is limited to low energies and masses, nevertheless, combined with dispersion relations and elastic 
unitarity it leads to a successful description of meson dynamics up to energies around 1 GeV, generating resonant states 
not originally present in the lagrangian, without any a priori assumption on their existence or nature. In particular, we 
find the p and a resonances as poles on the second Riemann sheet of nn elastic scattering amplitudes. With this 
approach we can then study some of these resonances properties, like their spectroscopic nature through their mass 
and width dependence on N c , or their dependence on the pion mass in order to connect with lattice studies. In the 
following sections we review this "unitarized ChPT" approach, named the Inverse Amplitude Method (IAM) [3, 4, 5], 
and then apply it to study the leading 1 /Nc behavior and the chiral extrapolation of the p and a mesons. 

The p and a resonances appear as poles on the second Riemann sheet of the (I, J) — (1,1) and (I, J) = (0,0) %% 
scattering partial waves of definite isospin, / and angular momentum J, respectively. Elastic unitarity implies for these 
partial waves, t(s), and physical values of s below inelastic thresholds, that 



where s is the Mandelstam variable and p is the center of mass momentum. Consequently, the imaginary part of the 

inverse amplitude is known exactly. However, ChPT amplitudes, being an expansion t ~ ?2 +?4 H , with = 0(p k ), 

can only satisfy Eq. (1) perturbative ly 




a(s), 



with a(s)=2p/y/s, 



(1) 



lmt 2 (s) = 0, lmt 4 (s) = o(s)t2(s) . . . 



(2) 



and cannot generate poles. Therefore the resonance region lies beyond the reach of standard ChPT. This region 
however, can be reached combining ChPT with dispersion theory through the IAM [3, 4, 5]. 



The analytic structure of the %% scattering amplitude t (s), consisting on a right cut extending from s,h = Am 2 n to oo, 
and a left cut from — °o to 0, allows to write a dispersion relation for the auxiliary function G(s) = t\ (s) /t{s) 

G(s)=G(0)+G'(0)s + iG"(0)s 2 + - / ds' > +LC(G)+PC, (3) 

where the integral over the left cut has been abbreviated as LC(G) and PC stands from possible pole contributions 
corresponding zeros of t . The different terms in Eq.(3) can be evaluated using unitarity and ChPT as follows: The right 
cut can be exactly evaluated taking into account the elastic unitarity conditions Eqs. (1), (2), ImG(i') = — a{s')t 1 {s') = 
— I1TU4 (s' ) , for s' G (4m^ , 00) . The subtraction constants only involve the amplitude and its derivatives evaluated at s = 0, 
so they can be safely approximated with ChPT: G(0) ~ r 2 (0) - f 4 (0), G'(0) ~ 4(0) - ^(0), G"(0) ~ -t'l (0). The left 
cut, which is suppresed by 1 /s' 3 (s' — s), is weighted at low energies, so it is appropriate to approximate it with ChPT: 
LC(G) ~ — LC(t4). The pole contribution only appears in the scalar wave, which vanishes at the so called Adler zero. It 
counts 0(p 6 ) and and has been calculated explicitly [6] and is not just formally suppressed, but numerically negligible 
except near the Adler zero, away from the physical region. 

Neglecting PC for the moment, and taking into account that t2 (s) is just a first order polynomial in s, and that a 
dispersion relation can be also written for t\, we can write Eq.(3) as 

G{s) = 1$ ~ f2(0) " ms ~ f4(0) " tm ' s W (0) " 2 " I £ ^ Jtf'-s-ie) LC(f4) = tl[s) U{S) > (4) 

which immediately leads to the IAM formula t IAM (s) = t2 u\^t 4 u\ ■ The IAM formula satisfies exact elastic unitarity and, 
when reexpanded at low energies, reproduces the ChPT expansion up to the order used to approximate the subtraction 
constants and the left cut. Here we have presented an 0(p 4 ) IAM but it can be generalized to higher chiral orders. 
Note that in the IAM derivation ChPT has been always used at low energies, to evaluate parts of a dispersion relation 
whose elastic unitarity cut has been taken into account exactly. Thus, there are no additional model dependencies in 
the approach, which is reliable up to energies where inelasticities become important. Taking the pole contribution into 
account leads to a modified IAM formula [6] which is almost indistinguishable from the ordinary one except in the 
Adler zero region, where it fixes some problems of the ordinary IAM with the Adler zero [6]. Actually, this modified 
IAM formula is the one used in this work since, as it will be shown below, one amplitude pole gets near the Adler zero 
region. 

This simple IAM formula is able to reproduce nit scattering phase shift data up to roughly 1 GeV and generates 
the poles associated to the p and a resonances with values of the LECs compatible with standard ChPT [5]. The 
l/N c expansion is implemented in ChPT through the LECs, whose leading l/N c scaling is known from QCD. Also, 
the quark mass dependence implemented in the IAM agrees with that of ChPT up to the order used. Hence, it is 
straightforward to study the leading 1 /Nc behavior and the m dependence of the resonances generated with the IAM, 
which we proceed to expose in the following sections. 

NATURE OF RESONANCES FROM THEIR LEADING 1 /N c BEHAVIOR 

The QCD 1 /N c expansion [2] provides a clear definition of qq bound states: their masses and widths scale as (9(1) and 
0(1 /N c ) respectively. The QCD leading 1 /N c behavior of the ChPT parameters (f„, m n and the LECs) is well known. 
Hence, by scaling with N c the ChPT parameters in the IAM, the N c dependence of the p and a mesons mass and width 
has been determined [7, 8]. They are defined from the pole positions as ^fs~^f e =M—iY. Note that we should not take 
too large Nc values, since the N c — )■ °° is a weakly interacting limit, where the IAM approach is less reliable [9]. Also, 
for very large N c values, even a tiny admixture of qq in the physical state would become dominant, but this does not 
give any information about the dominant component of the = 3 physical state. 

Fig. 1 (left) shows the p and a mass and width N c scaling. It can be clearly seen that the p follows remarkably well 
the expected behavior of a qq state, confirming that the method can obtain the correct N c behavior of well known qq 
states. In contrast,the a does not that qq pattern, allowing us to conclude that the a dominant component is not qq. 

Loop contributions play an important role in determining the a pole position. Since they are 1 /N c suppressed 
compared to tree level terms, it may happen that for larger N c values they become comparable to tree level 0(p 6 ) 
contributions, which are subdominant in the chiral expansion, but not N c suppressed. Thus we checked the 0(p 4 ) 
results with an 0(p 6 ) IAM calculation [8]. We defined a j 2 -like function to measure how close a resonance is from a 




FIGURE 1. Left: p and a 1 /N c scaling 0(p 4 ). Right: a l/N c scaling 0(p 6 ). 



qq behavior. First, we used it at 0(p A ) to show that it is not possible for the a to behave predominantly as a qq state 
describing simultaneously the data and the p qq behavior. Next, we obtained an 0(p 6 ) data fit where the p qq behavior 
was imposed. Figure 1 (right) shows the M a and Y a N c scaling obtained from that fit. Note that both M a and Y a grow 
nearNc = 3, confirming the 0(p A ) result of a non qq dominant component. However, for N c between 8 and 15, where 
we still trust the IAM, M a becomes constant and T a starts decreasing. This may hint to a subdominant qq component, 
arising as loops become suppressed as N c grows. Finally, by forcing the a to behave as a qq, we found that in the best 
case this subdominant component could become dominant around N c > 6 — 8, but always with an N c — » °° mass above 
1 GeV instead of its physical <~ 450 MeV value. This supports the emerging picture of two low energy scalar nonets, 
one of exotic nature below 1 GeV and another of ordinary qq nature above 1 GeV. 



CHIRAL EXTRAPOLATION OF THE p AND a RESONANCES 

ChPT also provides an expansion of m K in terms of m (at leading order m\ <~ m). Thus, by changing m n in the 
amplitudes we see how the IAM poles depend on m. We report here our analysis of the p and a properties dependence 
on m K [10]. 

The values of m K considered should fall within the ChPT applicability range and allow for some elastic regime 
below KK, that would almost disappear if m n > 500, which would be the most optimistic applicability range. We 
expect higher order corrections to be more relevant as m % increases. Thus, our results become less reliable as m n 
grows. 

Fig. 2 (top left) shows the evolution of the a and p pole positions as m K is increased. In order to see the pole 
movements relative to the nn threshold, which is also increasing, we use units of m K , so the threshold is fixed at 
^/s = 2. Both poles move closer to threshold and they approach the real axis. The p poles reach the real axis at the 
same time that they cross threshold. One of them jumps into the first sheet and becomes a bound state, while its 
conjugate partner remains on the second sheet practically at the very same position as that in the first. In contrast, the 
a poles go below threshold with a finite imaginary part before they meet in the real axis, still on the second sheet, 
becoming virtual states. As m n increases, one pole moves toward threshold and jumps through the branch point to the 
first sheet staying in the real axis below threshold, very close to it as m K keeps growing. The other a pole moves down 
in energies away from threshold and remains on the second sheet. Similar movements were found within quark models 
[12] and finite temperature and density analysis [13, 14]. 

Fig. 2 (top right) shows our results for the p mass dependence on m K compared with some lattice results [15], and 
the PDG value for the p mass. Now M p is defined as the point where the phase shift crosses 7T/2, except for those m K 
values where the p becomes a bound state, where it is defined from the pole position. In view of the incompatibilities 
between different lattice collaborations, we find a qualitative good agreement with lattice results. The M p dependence 
on m n agrees also with estimations for the two first coefficients of its chiral expansion [16]. 

In Fig. 2 (bottom left) we compare the m, dependence of M p and M a , normalized to their physical values. The 

bands cover the LECs uncertainties. Both masses grow with m n , but M a grows faster than M p . Above 2.4m^ hys , we 
show two bands since the two a poles lie on the real axis with two different masses. 

In the bottom center panel of Fig. 2 we compare the m K dependence of T p and T a normalized to their physical 
values: note that both widths become smaller. We compare this decrease with the expected phase space reduction as 




FIGURE 2. Top Left: Movement of the <J (dashed lines) and p (dotted lines) poles for increasing m n (direction indicated by 
the arrows) on the second sheet. The filled (open) boxes denote the pole positions for the a (p) at pion masses = 1,2, and 

3 x m^ s , respectively. Top Right: Comparison the IAM M p dependence on m % with some recent lattice results[15]. Bottom Left: 
Comparison of the p (light) and o" (dark) mass dependence on m n . Bottom Center: Comparison of the p (light) and cr (dark) width 
dependence on m n . The dotted (p) and dot-dashed (cr) lines show the decrease due to only phase space assuming a constant coupling 
to Tin. Bottom Right: p and a couplings calculated from the pole residue. In all panels, the bands cover the LECs uncertainty. 



resonances approach the nn threshold. We find that T p follows very well this expected behavior, which implies that 
the p nn coupling is almost m K independent. In contrast, r<j deviates from the phase space reduction expectation. This 
suggests a strong m K dependence of the a coupling to two pions, which we confirm with a explicit calculation of the 
resonances couplings from the pole residues as shown in the bottom left panel. 
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